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Abstract 

This note provides explicit techniques to compute the exponentials of a variety of anti- 
Hcrmitian matrices in dimension four. Many of these formulae can be written down directly 
from the entries of the matrix. Whenever any spectral calculations are required, these can be 
done in closed form. In many instances only 2x2 spectral calculations are required. These 
formulae cover a wide variety of applications. Conditions on the matrix which render it to 
admit one of three minimal polynomials are also given. Matrices with these minimal polyno- 
mials admit simple and tractable representations for their exponentials. One of these is the 
Eulcr-Rodrigues formula. The key technique is the relation between real 4x4 matrices and 
the quaternions. 

PACS Numbers: 03.65. Fd, 02.10. Yn, 02.10.Hh 



1 Introduction 

Finding the exponential of a 4 x 4 anti-Hcrmitian matrix explicitly is a problem which is of 
importance to quantum physics and its applications, especially to quantum optics, quantum 
information processing and computation. The problem of computing the solution to 

U = iH{t)U, U(0) =h,U£ U{A) 

with H(t) a 4 x 4 Hermitian matrix, arises in the study of four -level (or "two-qubit") systems. 
The solution to this problem can be reduced to the problem of computing the exponential 
exp(i~/(i)H) with H a 4 x 4 Hermitian matrix, typically different from H(t), and ^(t) some 
function. This reduction is achieved via either a passage to a rotating frame, approximations 
such as the rotating wave approximation, or techniques such as the Wei-Norman expansion 
or the Magnus expansion, or a combination thereof, [19, 21, 15]. Further, in the context of 
controlling four-level quantum systems, it is known that the unitary generators obtainable 
by allowing arbitrary time-varying external fields is precisely those obtainable by using only 
piecewise constant fields. This is the "controllability with admissible controls is equivalent 
to controllability with piecewise constant controls" result of [11], valid for any compact Lie 
group. Now determining the unitary generator after the application of a constant field to a 
four-level system is a matter of exponentiating 4x4 anti-Hermitian matrices. Further impetus 
to this question is provided by the issue of universality in quantum computation. Indeed, in 
quantum computation, due to various universality type results, it is known that to synthesize 
any quantum circuit it suffices to realize unitary matrices which are the tensor products of the 
identity matrix and unitary matrices which have either size two or four [see, for instance, [4]]. 
This, in turn, is equivalent to generating unitary matrices of size two or four selectively (this 
is essentially what the identity matrices in the tensor product amount to - the identity refers 
to the fact that the external interaction which seeks to address certain levels or qubits, does 
not disturb the other qubits, i.e., the addressing of the target qubits or levels is selective). In 
particular, this requires the accurate computation of exponential of anti-Hermitian matrices 
of size two or four. 

The purpose of this note, keeping the above goal in mind, is twofold. First we point out 
that the formulae in [18] extend in a straightforward manner to provide explicit closed-form 



formulae for the exponentials of a variety of matrices in su(4) - the Lie algebra of 4 x 4 anti- 
Hermitian matrices with null trace. These formulae already cover a wide variety of physical 
applications. Secondly we characterize when a matrix in su(4) admits either a quadratic 
minimal polynomial or a Euler-Rodrigues type formula for its exponential. In either instance 
the exponential of the matrix has a particularly simple representation. 

It is obvious that there is no loss of generality in assuming that the matrix being ex- 
ponentiated has zero trace. It is noted further that, in most instances, in the problem of 
exponentiating ^2^ =1 "f(t)H , H € su{4), one can assume 7 to be constant. To illustrate 
this consider the well-known formula, exp(ia(t)I 2 <8> cr x + ib(t)I 2 <S> <J y + ic(t)l2 <S> o z ) = 
cos(A(i))/ 4 + ^^-{ta(t)I 2 (g>a x + ib{t)I 2 ®a v + ic(t)I 2 <8> <r z ),X{t) = ^/ a{t) 2 + b(t) 2 + c(t) 2 . 
This formula would not suffer any modifications, beyond X(t) being autonomous, if each of 
a(t), b(t),c(t) were constant. In particular all the results of Section 3 extend verbatim to the 
case where the j(t) is not constant. 

The formulae provided in this note rely on an associative algebra isomorphism between 
H £g> H and gl(4,R), with H standing for the skew-field of the quaternions and gl(4,R) 
representing 4 x 4 real matrices. This isomorphism, known from the theory of Clifford Algebras 
[12], has only recently been used in concrete (numerical) linear algebra questions. To the best 
of our knowledge the innovative work of [9, 13, 6, 14] on the eigenvalue problem for a variety 
of structured real matrices is the first such instance. 

In [18] the same isomorphism was used to compute closed-form formulae for exponentials 
of structured real matrices. This is indeed, the point of departure for this note. Given the 
close relationship between the quaternions and the Pauli matrices, it seems plausible that 
the basis for gl(4,R) provided by the aforementioned associative algebra isomporphism is 
essentially the basis for u{4) provided by the various Kronecker products of the Pauli matrices 
and I 2 . It is tempting to believe that this correspondence is as elementary as assigning, for 
instance, the elements i(g)l,j®l,fc(g)l in H ® H to the matrices a x ® I 2 , a y ® I 2 , a z ® 1 
etc., However, a moment's reflection shows it cannot be this simple. In light of this, it is 
a pleasant circumstance that the aforementioned plausible connection is indeed valid. The 
precise correspondence is presented in the next section. 

It is worth mentioning that the results presented here can also be used to exponentiate 



matrices in su(3). Indeed, one has to just embed such a matrix as a principal submatrix 
in a 4 x 4 matrix, with the rest of the 4x4 matrix consisting of zero entries. A different 
application would be to exponentiate matrices in so(6, R) (the Lie algebra of 6 x 6 real anti- 
symmetric matrices) matrices. S'f/(4) is a double cover of SO(6,R), [12]. By making this 
covering homomoprhism explicit, one can reduce the problem of exponentiation in so(6, R) 
to finding exponentials in su(4). 

The balance of this note is organized as follows. In the next section, the relation between 
the H ® H basis for gl{&, R) and the Pauli tensor product basis is made explicit. The same 
section also establishes notation used throughout this note. The third section presents su(A) 
analogues of the results of [18]. In particular, several illustrations drawn from important 
applications are given. The fourth section presents conditions which ensure that a sw(4) 
matrix either has a quadratic minimal polynomial or admits a Eulcr- Rodrigues' formula. 
The same section also presents conditions equivalent to a su(4) matrix to stem from a normal 
matrix (i.e., if X = B + iC, then BC = CB, where B,C are the skew-symmetric and 
symmetric parts of A e su(4), respectively). The final section offers conclusions. 

2 Notation and Preliminary Observations 

The following definitions and notations will be frequently met in this work: 

• gl(n, R) represents the algebra of real nxn matrices. This is, of course, the Lie algebra 
of the Lie groups of real invertible matrices. 

• SU(n) represents the Lie group ofnxn unitary matrices of determinant one. su(n) 
represents the corresponding Lie algebra ofnxn skcw-Hermitian, traceless matrices. 
su(2) (g> su(2) is the Lie algebra spanned by matrices of the form 7 2 <£> U + V <X> h , U, V G 
su(2). Note it is customary to use the terminology "anti-Hermitian" for skew-Hermitian 
matrices. 

• R n represents the matrix with 1 on the anti-diagonal and zeroes elsewhere. Any nxn 
matrix A, satisfying A T R n + R n A — said to be per skew symmetric. Per symmetric ma- 
trices are those matrices, A, which satisfy X T R n = R n X. Such matrices are symmetric 
about the anti-diagonal. 



/ 

• J2n is the 2n x 2n matrix which, in block form, is given by J 2n = 



y In On J 

Matrices, Z, satisfying Z T Jm — JinZ are called skew-Hamiltonian (sometimes anti- 
Hamiltonian) . The term "Hamiltonian" will not typically be used in the sense of quan- 
tum mechanics, unless specified to the contrary (i.e., it will not be used to mean a 
Hermitian matrix). 

• Throughout H will be denote the field (more precisely the division algebra) of the 
quaternions, while P stands for the purely imaginary quaternions, tacitly identified 
with R 3 . Further, in this note the symbol i will be used for both the corresponding 
complex number and the corresponding quaternion. The context should make it clear 
which of the two is implied. 

Remark 2.1 • i) Throughout this note, use of the following observation will be made: 
Let X be an n x n matrix satsifying X 2 + c 2 I n = 0, c ^ 0. Then e x = cos(c)/„ + ^lx. 
Here c 2 is allowed to be complex, and c is then taken to be \pre 1 ^- , with c 2 = re 10 ,9 e 
[0,2tt). 

• The fact that any matrix which satisfies X 3 = —c 2 X, c ^ 0, satisfies e x = I + Sln ^ X + 
i-cos(c) ^2 ^ ne Euler - Rodrigues's formula) will also be used. Once again c 2 is per- 
mitted to be complex. Note that any matrix which satisfies X 2 + c 2 I n = 0,c ^ 0, 
automatically satisfies X 3 = —c 2 X, c ^ 0. For such matrices the exponential formula 
in i) is better to work with than the Rodrigues' formula. Therefore, in this note we will 
allude to a matrix admitting an Euler- Rodrigues formula only if its minimal polynomial 
is of the form x 3 + c 2 x. 

• in) Explicit formulae for e A can be produced if the minimal polynomial of A is known 
and it is low in degree (cf., [1] where such formulae are written down from the char- 
acteristic polynomial). However, since the corresponding explicit formulae for e A are 
more complicated than the ones corresponding to i) and ii), they will not be pursued 
here. 



H®H and gl{&, R): The algebra isomorphism between H ®H and gl(4, R), which is central 
to this work is the following: 

• Associate to each product tensor p (g> q £ H ® H , the matrix, M p ® q , of the map which 
sends x £ H to pxq, identifying R 4 with H via the basis k}. Thus, if p = 
Pa + Pli + Pij + P3k; q = q + q 1 i + q 2 j + q a k, then 

M pm = [x\y\u\v] 

with x, y, u, v, the columns of the matrix M p ® q , given by the vectors in R 4 representing 
the quaternions pq,piq,pjq,pkq respectively. Here, q — qa — q\i — q 2 j — q^k. 

• Extend this to the full tensor product by linearity, This yields an algebra isomorphism 
between H®H and 5/(4, R). In particular, a basis for gl(4, R) is provided by the sixteen 
matrices M Ci8e as e x , Cy run through 1, i 1 j, k. 

This connection, which is known from the theory of Clifford Algebras, has been put to 
great practical use in solving eigenvalue problems for structured matrices by Mackey et al., 
[9, 13, 6, 14]. It can also be used for finding exponentials, e A ,A £ gl(4,R), [18]. 

Remark 2.2 Canonical Form for X £ su(4) : Let X = iH, with H Hermitian and traceless. 
Then 

3 3 33 

H = ^ Oiih <E)0-i+^2 Pi a i ® J 2 + X! 7 -? fe °J ® U k,Oii,(3i, Jjk € R (2.1) 
i— 1 i— 1 j — 1 fc— 1 

It is well known that via conjugation by a local unitary transformation (ie., conjugation via 
&U £ SU(2) <g> SU{2)) H can be put into the form 

333 

di h ® o-j + ^o-j j 2 + m (2.2) 

i— 1 i— 1 i— 1 

with di,bi,Ci £ R. We will use this canonical form at some points in Section 4 (but not 
in Section 3). Furthermore this local unitary transformation is determined by finding the 
singular value factorization of the real 3x3 matrix (jjk)- But this amounts to finding the 
spectral factorization of a real 3x3 symmetric matrix - which can be performed in closed 
form, [3]. 



Relation to the Pauli Tensor Product Basis As mentioned in the introduction, the 

above basis for gl{A, R) is closely related to the basis Ui ® aj, i, j = 1, . . . , 4(with cr = ^2,0"! = 

o~ x , &2 — &y, 03 = CTfe). The precise relation is tabulated below: 
Pauli Tensor Basis Quaternion Tensor Basis 
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3 Some Closed Form Formulae for Exponentials in sw(4) 

In this section, we provide closed form formulae for the exponentials of several matrices in 
su(4), without resorting to the canonical form in Equation (2.2). These formulae are based on 
expressing the matrix in question as a sum of commuting summands, each of which satisfies 
the condition in i) of Remark (2.1). These formulae can be divided into two classes: i) those 
which can be directly written down from the entries of the matrix; ii) those that require 
the spectral factorization of an associated real 3x3 symmetric matrix. This latter spectral 
factorization can be achieved in closed form, [3]. In particular, for several cases only a 2 x 2 



spectral factorization is needed. These will be pointed out. Since most of these formuale are 
the su(4) analogues of the results in [18], proofs will be provided only for cases not considered 
in [18]. In the interests of brevity, we have not considered analogues of every possible result 
in [18]. 

Remark 3.1 Consider X e sw(4), written as X = B + iC, with B,C real. Suppose it is 
skew-Hamiltonian, for instance. Then a simple calculation reveals that the real matrices B, C 
are skew-Hamiltonian as well. Hence so is the real matrix B + C. This observation yields 
the H ® i? representation of such anl€ su{A). The basic properties used in exponentiating 
the corresponding real matrix B + C in [18] was that it could be expressed as the sum of 
commuting summands, each of which is annihilated by a polynomial of the type in i) of 
Remark (2.1. Now these properties are not vitiated by the presence of the imaginary unit 
i m X. Therefore their exponentials are similarly found. The only difference is that the 
hyperbolic trigonometric functions in the formula for e B+c will now be replaced by their 
ordinary trigonometric equivalents. Similar arguments hold if X is perskewsymmetric etc., 

3.1 Exponentials Directly From the Entries 

Below a list of three families of matrices in sw(4), whose exponentials can be directly found 
from their H ® H representations, is presented. 

1. Symmetric, Tridiagonal, Su = Consider 

/ \ 

a O 



S = i 



V 
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7 

Since such matrices arise in several applications, it is interesting to note that they can 
be easily exponentiated. Indeed, such an S has the following representation 

S = i[M pS)l + M q ®j + M rm ] =X + Y + Z 

with p = (0, |, 0), q = (|,0, ^±2), r = (0, ^,0), a, 0, 7 e R. In terms of the Pauli 



tensor basis, S is i§((T x ® a x ) +i^a y <8>cr y -M^^^ ®cr x +i^^-a z <&a x . Now note that 



Y commutes with both X and Z, while X and Z anticommutc. Further each squares 
to a negative constant times the identity. So e s is given by 

e b = [cos(Ai)/ 4 + i , (M p8i + M r8fc )][cos(A 2 )/ 4 + « — r-^Mj®j] 
Ai A 2 

In terms of the Pauli matrices this becomes e s — [cos(Ai)/ 4 + i ™^ 1 ^ (^(a x ® a x ) + 

^(T z m x )][cos(X2h+i 9j ^ 1 (^ y ^a y +^I 2 m x )} withAj = § n//? 2 + (7 - ") 2 , A 2 = 

V/5 2 + (7 + ") 2 - 

2. Perskewsymmetric X : Such an X e su(4) satisfies, in addition, Jf T i? = —RX. 
Such matrices are expressible in the form 

i[pia z <g> h + P2&x <8> cr z + oea y <g> a z + q\I 2 ®o z + <72Cz ® &x + P&z <8> cr y ] 
Their exponential is given by 

r /x w ,sin(Ai. T Mr /-x nt .sin(A 2 . T 

[cos(Ai)i4+« — : (pit7 z (8)/ 2 +p 2 tT a; (8)CT z +a(T y (g)cr z )J[cos(A2)i4+« — r {qih®cF z +q2(T z ®a x +pa z ®ay)\ 

Ai A 2 

with A x = v/|M| 2 +a 2 , A 2 = Vll9ll 2 +/3 2 

3. Skew-Hamiltonian X: These matrices satisfy, in addition, X T J = JX. Such matrices 
are associated with time-reversal symmetries, [8]. More specifically, a Hamiltonian (in 
the usage of quantum mechanics), H, i.e., a Hermitian H, is associated to time- reversal 
symmetry if H T J = JH. Clearly if H satisfies this additional condition, so does 
X = iH. Such matrices are expressible in the form 

i[bh + Picry <E> o y + p 2 / 2 ® a z + p^h ®o x + ca z <g> a y + da x ® a y \ 

Their exponential is given by 

e lj) [cOs(A)/4+Z bm ^ {pi(Jy®(Jy+P2l2®<?z+Psh®Gx+CCrz®Gy+d<J x ®<Jyj\, A = VII P II 2 + c2 + rf2 
A 

3.2 The Purely Imaginary Case 

The following algorithm for exponentiating a matrix X e su(4), which is simultaneously 
symmetric (equivalently purely imaginary) follows directly from the corresponding algorithm 
for exponentiating purely real symmetric matrices in [18]. The only difference is that the 
coshQ, sinhQ in [18] will be replaced by cos(),sin(). Note that such an S will not have any 



terms in Equation (2.1) corresponding to members of the Pauli tensor basis, which contain 
precisely one a y term. 

• Represent the given symmetric S € su(4) as the matrix as i[Mp^i+M q 0j+M r k ],p, q, r e 
P. 

• Identifying the pure quaternions p, q, r with vectors in R 3 , find the spectral factorization 
of the real 3x3 symmetric matrix X T X, where X = [p \ q \ r]. Thus X T X = 

E3 9 T 

• Compute Ui = Xvi (Note Ui are almost the left singular vectors. The only difference is 
|| m ||= Gi). Then S = «X)Li M u z ® Vi - Hence, 

e s = H3 =1 (cos(a,)/ 4 + i^^-M^) (3.3) 

Definition 3.1 Bisymmetric Type: For several important examples only a 2 x 2 spectral 
factorization is needed (which is extremely easy to write). Since the archtypical example 
is provided by a matrix in su(4) which is, in addition, bisymmetric (i.e., simultaneously 
symmetric and persymmetric), we will, to avoid circumlocution, call all such matrices of the 
bisymmetric type. 

3.3 Illustrative Examples 

We provide some important illustrations of the formulae developed in this section. 
Illustration 1: Rabi Oscillations in Four Level Systems In [7] a detailed calculation, via 
a calculation of eigenvectors and eigenvalues, is provided to calculate the evolution of a four 
level system, being irradiated by three laser fields, under the rotating wave approximation and 
under resonance. Specifically, they consider a four level system with energy levels {Ek, k = 
1, . . . , 3} which satisfy E\ — E > E 2 — E\ > ... > E s — E 2 . This system is irradiated by three 
laser fields with frequencies uj k = E k — E k _i, k = 1, . . . , 3. After passage to a rotating frame, 
and under the assumptions of resonance and the rotating wave approximation, the unitary 
generator in the rotating frame satisfies 

iU = (E I 4 + C)U (3.4) 



with 
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Here the gi are the amplitudes of the three laser fields. Thus U(t) = e~ lEat exp{—iCt). In 
[7] exp(-iCt) is calculated by a direct calculation of the eigenvalues and eigenvectors of the 
matrix —iC. Now, — iC is precisely a symmetric, tridiagonal matrix with a zero diagonal - 
i.e., of the type considered in item 5) of the list. Thus, exp(-iCt) may be found directly and 
is equal to 



-iCt 



== [cos(Ai)Zi+i 



.sin(Ai) 



sin(A 2 



■(f3a y ®a y +('y+a)l2®c 



i A 2 
with a = -± gi t,f3= -\g 2 t, 7 - -^ 3 t, Ai = V/? 2 + ("-7) 2 , A 2 = ^/J 2 + ( 7 + a) 2 . 

A laborious but straightforward calculation confirms that the matrix entries provided by 
the above representation of exp(-iCt) coincide with those in [7]. 
Illustration 2: Josephson Junction In [22, 20] the following system is considered; 
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In [22] E o, E\o,Eji,Ej2 are taken to be constants reflecting current technology. Thus U(t) = 
e -iHt_ Now note that 

-iH = -i[^{E 00 + E w )h - ~E J2 CT X ®I 2 - ^Ejxh <g> u x + ^(E 00 - E w )a z <g> a z ] (3.5) 

In terms of the Pauli tensor basis this is —i[^(E 00 + E w )M im - \Ej 2 Mi® k - \Ej\Mk®j + 
\{E m - E 10 )M mk ] Hence, e -* Ht = e -i(E aa +E w )t e -im ^ with = -i[M pm + M q<S)j ], with 
the purely imaginary quaternions of the form p = pii+p 3 k, q = q 3 k. Thus, the singular value 
factorization of the 2x2 matrix 

<» ^ 

P3 <73 



has to be found. Thus, this is an example of the bisymmetric type. Specifically, the calculations 
proceed as follows: H = — i[M Ul ® Vl + M U2 ® V2 ) with v\ — cos 8i — sin6*fc, v 2 — sindi + cos 8k, 
Here tan(2#) = qV q _pT p - Further u\ = pi cos 8i — (p 3 + q 3 )s'm8k,u 2 — pis'mdi + (p 3 + 
qs) cos 8k. Then 1 1 v% | |= 1, i = 1, 2, while 1 1 m | |= ai = \J q T q cos 2 8 + r T r sin 2 8 — q T r sin(2#) , 



1 1 u 2 1 1 = 02 = Y Q T <7 sin + r T r cos 2 + g T r sin(2#) . Hence, 

. sincri,, lr sin<72,, , ,„ „. 

e lHt = [coso-i/-z M u 1 ® v J [cos a 2 I - 1 M U2 ® V2 \ (3.6) 

(Ti (72 

This reads, in terms of the Pauli matrices, as the following 

Illustration 3. Scalar Coupling Hamiltonian The matrix being exponentiated is X = 
—iH = i[ali + ba z <g> I 2 + cl 2 ® ct z + c?ct z ® cr z + e<r x ® a x + fa y <g> ct^]. This is the so-called 
scalar coupling Hamiltonian, and is widely used in NMR spectroscopy. The corresponding 
H <£) H representation is given by 

X = -i[aMi® 1 +bM m +cM j ®j+dM k ® k +eM j ® i +fMi ls , j ] = -i[aMi m +M p ®i+M q ®j+M r ® k ] 

with p = (b, e, 0), q = {f, c, 0), r = (0, 0, d). 

Hence e tx = e lat exp — it{M p ®i + M q ®j + M r ® k ) . Thus, it remains to find the exponential 
of the purely imaginary symmetric matrix —it(M p ®i + M q ®j + M r ® k ). Now notice that the 
real matrix — t[p \ q \ r], is a block diagonal matrix, with the northwest block a real 2x2 
matrix and the southeast block, the lxl matrix (—td). Hence, one needs to only find the 



singular value factorization of the real 2x2 matrix 

(" ') 

-t =]p\ Q] 

[e c) 

Hence this is also of the bisymmetric type. The corresponding right singular vectors of — t\p 



q\r], written as quaternions, arevi = cos8i—sm8j,v 2 = sm8i+cos8j,Vs = k, Heretan(2#) = 
qTg-pTp = fi^-^-e^ ■ Furtner "i = -t(bcos8i - / sin 8j), u 2 = -t(esin8i + ccos6j),u 3 = 
—tdk. Then 1 1 Vi \ \= 1, i = 1,2,3, while 1 1 m 1 1= 01 = t^J p T p cos 2 8 + q T q sin 2 8 — p T q sin(28), 
1 1 u 2 1 1= <T2 = t\jp T p sin 2 8 + q T q cos 2 8 +p T qsm(28), \ \ u 3 \ \= td Hence, 

CTl (T 2 (T 3 

Remark 3.2 There are several other practical applications which lead to the problem of 
exponentiating su(4) matrices of the bisymmetric type. Examples include superconduct- 
ing circuits for solid-state quantum computation [5], J cross polarization experiments [17], 



Heiscnbcrg Hamiltonians (under the assumption that only one of the three components of the 
magnetic field, assumed to be constant in time, is active during any period of time). 

Remark 3.3 The number of matrices which can be easily exponentiated in this fashion can be 
expanded by combining the above observations together with some useful conjugations. Two 
classes of such conjugations immediately spring to mind. The first is obviously the class of local 
unitary transformations. Thus, for instance the matrix X\ = i(aa z <£> a z + ba y ®h + cIi <g> <r y ) 
is explicitly locally unitarily equivalent to X2 = i{aa z <g> a z + ba x ® I2 + cl<i ®a x ). The former 
is not a symmetric matrix, while the latter is (in fact, it is of the bisymmetric type). However, 
e Xl is easily found once e X2 is. The second type of conjugation is via the so-called magic 
basis matrix (see [16], for instance). Explicilty, letting 
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it is known that V{so(A, R))V* = su(2) ® su{2). It is instructive to examine its effect on some 
of the other matrices considered here. Thus, for instance 
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/ 



• X, symmetric, tridiagonal with Xu = implies VXV* = i[^a y 0/2+ 2 a v ® a x + 
- L ^ L cFy ® a z + ^j 1 cr z <S) <J y ] If one writes such a matrix explicitly, it is not clear that 
it too can be written as the sum of commuting summands, each of which are easily 
exponentiated. 

• X, skew-Hamiltonian implies VXV* = i[al 2 ® h + ba y <£> a x + ca y ® a y + da y ® a z + 

ea x (g) I2 + jo z ® I 2 ] 

• X, perskewsymmetric implies VXV* = i[aa x ® a x + ba x ® a z + c/2 ® <J y + da y ® a y + 

ea z ®a y + ja x ® I 2 ] 



Thus, all such matrices are readily exponentiated. 



4 Euler-Rodrigues Formulae 

In this section we will find conditions which imply that a given X e su(4) admits one of 
three minimal polynomials. For such X the corresponding formulae for e x is very handy. In 
particular, one of them is an Euler -Rodrigues type formula, which explains the title of the 
section. Finally, we will provide conditions that an su(4) matrix X = B + iC is of the normal 
type, i.e., BC = CB. Note that, since X is a complex matrix, being of the normal type is not 
the same as being normal. For some of the results below we will appeal to the canonical form 
for X in Equation (2.2). This section will make use of the eigenvalue structure of matrices in 
su(4). However, one does not need to determine the eigenvalues themselves. 

Given an annihilating polynomial p for any matrix X , one can use p to find e x . There 
are at least two manners in which to achieve this. One finds the zeroes of p (but not the 
eigenvectors of X), and then proceeds to use any of a variety of methods (e.g, interpolation) 
to find e x [10]. Alternatively one can express, via p, higher order powers of X in terms of 
lower orders and use this to establish formulae for e x . In general, the representations of 
e x obtained via either method are not always easy to work with. There are however three 
instances when either method produces the same representation and this is particularly easy 
to manipulate. Specifically, these are: 

• 1. Quadratic Type I: p(x) = x 2 + c 2 . The corresponding formula for e x is e x = 
cos(c)/„ + ™&X. 

• 2. Quadratic Type II: p(x) = x 2 + 2(3x + j, f3 ^ 0. Now e x = e~ [(cos o + 13 si ^ ( a) )I + 
^X], with a = yj[3 2 - j. 

• 3. Cubic Type I: p(x) = x 3 + c 2 x. In this case, e x = I + s Me1x + l^£l X 2 

Remark 4.1 The results to be presented here should be seen as a complement to those in 
the previous section. In Section 3, only the first of the above minimal polynomials was used. 
There will be several X e su(A) which are amenable to the techniques of either section. 
Consider, for instance, X = —iJ(t)(a x a x + a y a y + a z a z ). This matrix arises in the study of 
quantum dots, [5]. The corresponding p is x 2 — 2iJx + 3J 2 . However, X is also the sum of 
three commuting terms, each of which is annihilated by a polynomial of the first type in the 
above list. 



We begin with an explicit expression for the characteristic polynomial of an X e su(4) in 
canonical form. 

Proposition 4.1 Consider anle su(4) in canonical form as given in Equation (2.2). Let 
its characteristic polynomial be x 4 + [ix 2 + vx + n. Then 

• i)M = 2ELi(«?+&? + c?) 

• ii) v = -8i(£? =1 a-ihci - n? =1 c;). 

• - = |{2(Eli(«f + &? + 4)? 4[(Ef =1 (o? + 6? + c 2 ) 2 

+ ^Eij^l^^j C 2 C 2 — ^Ey,tl:./;# dibiCjCk]} 

Proof: These formulae follow from Newton's identities, which imply that the coefficients 
of the characteristic polynomial can be expressed in terms of the trace of suitable powers of 
X, in conjunction with Tr (X) = 0. Further, Tr (X l ),i = 2, . . . , 4 were calculated by using 
the H ® H representation of X and looking for the 1 <£> 1 term in X 1 . It is worth emphasizing 
that the ease of quaternion multiplication renders it unnnecessary to calculate X 3 or X 4 fully. 
Indeed, besides calculating the 1 ® 1 term in X 3 , one needs to find only those terms in X 3 
which would yield al®l term in X 4 (and quaternion multiplication facilitates this process). 

We can now give a simple characterization of when X's minimal polynomial is of either 
quadratic type I or cubic type I. 

Proposition 4.2 X € su(4) has i) minimal polynomial p(x) — x 2 + c 2 iff v = and [i 2 = An; 
ii) minimal polynomial p(x) = x 3 + c 2 x iff v = = ir. Furthermore, in case i) c 2 = ^, while 
in case ii) c 2 = ji. 

Proof: First, in view of X's diagonalizability, p(x) — x 2 + c 2 is the minimal polynomial 
iff the characteristic polynomial has two distinct roots (which add up to zero) each repeated 
twice. Similarly, p(x) = x 3 + c 2 x is the minimal polynomial iff the characteristic polynomial 
has two simple distinct roots (which add up to zero) and a double root equal to zero. 

Suppose first that v = 0. Then the characteristic polynomial is a quadratic for x 2 . The 
first case occurs precisely when this quadratic has a double root, i.e., when fi 2 = An. Similarly, 



the second case occurs when one of the roots of this quadratic is nil, i.e., precisely when 7r = 
in addition. 

Conversely, suppose the minimal polynomial is p(x) = x 2 + c 2 . Now using the characteriza- 
tion of the coefficients of the characteristic polynomial in terms of the elementary symmetric 
functions of the eigenvalues, it follows that v = and [i 2 = An. Similarly if p(x) = x 3 + c 2 x 
the same characterization yields v = = n. 

Remark 4.2 i) Using these conditions it is easy to write down examples of X e su(A) 
which admit genuine Euler-Rodrigues formulae, i.e., X which have cubic Type I minimal 
polynomials. For instance, X — ® a x + a x ® I2 + <J y ® <J y + ca z (g> a z ), where c is any real 
solution of the quartic c 4 + 14c 2 — 8c+ 17 = 0. This quartic admits at least two real solutions. 
Indeed, if all solutions were complex, then they must be of the form a+ib, a—ib, —a+id, —a— id, 
since there is no c 3 term. It is easy to see that, if this is the case, then the coefficient of c 2 
has to be necessarily negative. Note further, that c, in this example, could easily be allowed 
to be time- varying. 

ii) It is noted in passing that one can write down the exponential of generic X which satisfy 
v = (i.e., those cases for which neither of [i 2 = A-k nor tt = hold), since in this case the all 
eigenvalues of X are distinct and the corresponding interpolation based formula [10] assumes 
a simple form. 

Characterizing when X has a minimal polynomial of Quadratic Type II via coefficients 
of the characteristic polynomials does not seem fruitful. Therefore, we provide a different 
characterization. For this characterization we do not require that X be placed in the form 
of Equation (2.2), though obviously the stated conditions would simplify for X in canonical 
form. 

Proposition 4.3 Let X e su(4) be expressed as M p ®\+M\® q +i[M r ®i+M s ®j+M t ®u\, with 
p, . . . ,t purely imaginary quaternions. Denote by C = [r \ s | t] . Then X admits x 2 + 2(3x + 7, 
with [3 ^ 0, as its minimum polynomial iff there is a f3 € R satisfying the following conditions: 

C T p = j3q (4.7) 
Cq - h 
Pq T - Go (C) = j3C 



where Co (C) is the matrix of cofactors of C. If these conditions hold, then i) 7 = — [|| p \\ 2 
+ ||g|| 2 + ||r|| 2 + ||s|| 2 + ||i|| 2 ];ii)/3 = ^. 

Proof: From a variety of viewpoints it should be clear that if x 2 + 2(ix + 7 is to be the 
minimum polynomial of X, then necessarily f3 is purely imaginary, while 7 is real. Using this 
fact, the above conditions stem from a direct calculation of X 2 . 

Remark 4.3 The purpose of this remark is to identify some situations, under which, the 
system of equations in Equation (4.7) admits solutions. 

i) When C has rank one, Equation (4.7) always possesses solutions, i.e., it is always possible 
to find p, q, $ satisfying them for the given C. Such a C always possesses a representation 
of the form C = uv T , with u T u = v T v (it is easy to find such a representation). Picking 
(3 = u T u = v T v,p = u^J~fi,q — v\Jfi we find that Equation (4.7) is satisfied. This yields 
a systematic procedure to construct examples admitting a quadratic minimal polynomial of 
type II. 

ii) Conversely, starting with a non-zero p one can find q, C such that Equation (4.7) always 
holds. The key to this is to beserve that if C is invertible in Proposition (4.3), then the 
conditions given in Equation (4.7) can be written in a different form, to wit: CC T p = f3 2 p, q = 
C~ 1 (f3p); f3(pp T — CC T ) = det(C)I. This yields a method to construct more examples of X 
admitting a quadratic minimal polynomial of Type II. Pick ap^O. Choose (3 = a/1 + p T p, 
and pick a C satisfying dei(C) = —/3 and CC T — I + pp T . Finally, set q = C~ 1 (f3p). 

This can always be achieved by picking C to be a solution of the equation CC T — I + pp T 
with a negative determinant, since I+pp T is positive definite and thus possesses a square root. 
For instance, one could multiply the easily determined unique positive definite square root of 
I + pp T by diag(l, 1, —1) to obtain a C with determinant — (3. Further, det(CC T ) = 1 +p T p 
and obviously CC T p — f3 2 p. Indeed, the eigenvectors of / + pp T are p (with eigenvalue 
1 + p T p = [3 2 ), and any two vectors orthogonal to p (corresponding to eigenvalue 1 with 
double multiplicity). 

iii) If precisely one of p or q is zero, then there is no solution to Equation (4.7). When 
both arc zero, there is a solution iff CC T is proportional to the identity matrix, i.e., iff 
the vectors r, s and t are orthogonal and have the same length. Note, in this case X = 



i\M\ + M 2 + M 3 ], with the Mj commuting, and each with a quadratic minimal polynomial of 
type I. Further, this is precisely the case wherein the canonical form X, as in Equation (2.2), 
is X = i\c\o x ® a x + c 2 a y ®a y + c^a z ® a z , with either ci = c 2 = c 3 (in the event det(C) > 0) 
or ci = c 2 = — c 3 ( in the event det(C) < 0). However, one does not require passage to this 
canonical form for finding e x . 

iv) Similar conditions can be written down one p,q,r,s,t for X to admit other minimal 
polynomials. We omit them in the interests of brevity. 

Conditions for "Normality" Next, given X = B + iC, we characterize, when [B, C] = 0, 
i.e., when the real matrix B + C is normal. The motivation should be obvious - it is possible 
to exponentiate both B and iC in closed form, and hence X. While the statement of this 
result uses the canonical form given by Equation (2.2), much of the proof does not require it. 

Proposition 4.4 Let X = B+iC = M p ®i + i® q +iM r ® i+s ®j +t ®k,p, . . . , t e P be in canonical 
form. Suppose, without loss of generality, that at least one of p, q is non-zero. Then [B,C] =0 
iff the following conditions hold: 

• i) p ^ 0, q = 0: ai = a 3 = c\ = c 3 = 

. ii)p^0,g^0: ai =a 3 = 61 =63 = 0,| £|=| £1=1 

• iii) p = 0, q ^ 0: b 1 = b 3 = c x = c 3 = 

Proof: For any X e su(4) (even those not in canonical form) a quick calculation reveals that 
[B, C] = 2[p x r <E i +p x s ® j +p x t <g> k + r ® (g x i) + s® (g x j) +t ® (g x k)]. If X is in 
the canonical form in Equation (2.2) then we have 

p = (-a 2 ,0,0);g = (0,6 2 ,0);r = (6 3 ,ci,0);s = (c 2 , a 3 , m), i = (&i,0,c 3 ) 

Hence [£>, C] is twice the matrix representation of b\b 2 i <£> i + (c 3 6 2 — cia 2 )fc ® i + aia2j <E 
j — a 2 a 3 fe (E> j — & 3 6 2 « <E fc + (a 2 c 3 — Ci6 2 )j <E fc. The conclusion follows from this. 

Remark 4.4 By applying the proof of the previous result to X not in canonical form, one can 
deduce other commutativity results. For instance, suppose Y is in canonical form, and one 
defines Y x = iQ^Li <E cr, + X)»=i ® ^2) and Y~ 2 = i X)»=i Ci^i ® Then [Yi, Y 2 ] = 
iff i) S3. = ^; ii) a = ^; and iii) ^ = T o sec this let X = B + iC = V*YV, with V the 



magic basis matrix [see Remark (3.3)]. Then [ii, I2] = iff B + C is normal. Note, while Y 
is in canonical form, X is not. 

5 Conclusions 

In this note, closed form formulae are provided for exponentials of several important anti- 
Hcrmitian x4 matrices. These matrices cover many important applications. The basic tech- 
nique is the isomorphism between real 4x4 matrices and H <g> H. We believe that this 
connection is aptly suited to exploit the properties of su(4) stemming from its direct sum 
decomposition into the real skew-symmetric matrices and the purely imaginary symmetric 
matrices. While no claim to the superiority of the representation of the exponential provided 
by this work is made, it is our hope that further research will yield more applications of these 
formulae. 
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